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ABSTRACT 



The work done in this investigation is part of a continuing program 
designed to develop probabilistic strength and life models for filament 
composite materials. In this particular experiment, groups of single fibers 
from two chemically identical different production spools of graphite were 
loaded to identical sustained tensile loads in order to determine their life 
statistics while under stress. While previous work has produced models 
of composite reliability as a function of strength, this work develops a 
model which will be necessary to predict the reliability of composite fibers 
as a function of time. Parametric and non-parametric methods were used 
to quantify the relationship between fiber reliability as a function of stress 
and as a function of time. 
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I. INTRODUCTION 



While still in its infant stages, the use of composite materials in 
manufacturing is steadily increasing, particularly in high reliability, low 
accessibility environments such as aircraft fuselages and satellite fuel tanks. 
High strength to weight ratios, combined with minimal maintenance 
requirements and extreme resistance to corrosion make composites ideally 
suited for structures previously built exclusively of steel or aluminum. 
However, because composite usage has such a short history, there is no 
significant data base which would allow accurate reliability predictions. 

Through an ongoing program of research and experimentation, The 
Advanced Composites Laboratory at the Naval Postgraduate School at 
Monterey, CA, is defining a probabilistic model which is based on the 
failure mechanism of composites in tension. Small scale designs which 
exploit fiber strength under tension, have accurately verified the Local 
Load Sharing model. More recent work is concentrating on the application 
of this model to large scale structures. With this in mind, it is important to 
understand that a structure must not only be predictable (reliable) in 
strength, but also in life. Because of their non-homogeneity, the study of 
composite materials must be further broken down to the study of not only 
individual constituents, but also of their inter-relationships with each other. 
For our purposes here, all work is done exclusively with graphite/ epoxy 
composite.; therefore, a proper understanding of the fiber, the matrix, and 
the interaction between them , is essential to formulating an accurate model 
of composite reliability. 



The primary objective of this project was to examine a large sample 
space of graphite fibers under different load conditions, and formulate a 
model of their reliability in life. The fiber life results would then be 
applied to composite life experiments now on-going designed to predict the 
life of composite structures. With data from these two projects in hand, 
one could then ostensibly isolate the solution to the fiber/ matrix interface 
problem and produce a realistic life prediction model. 
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II. BACKGROUND 



A. STRENGTH DETERMINATION OF FIBER REINFORCED 
COMPOSITES 

Fiber Reinforced Composite materials generally consist of high 
strength fibers bonded together in a ductile matrix. In these materials the 
fibers carry essentially all of the stresses imposed on the structure while 
the matrix serves as the primary vehicle for transferring load between 
adjacent fibers, known as the Local Load Sharing model. In this model, 
the local load carried by a bundle of fibers is carried equally until the 
weak fibers begin to rupture. When a fiber breaks, its load is transferred 
to the adjacent fibers through the matrix. Rosen has modeled the load 
transfer by the existence of an ineffective length at the end of each 
(failed) fiber [Ref. 1]. In this region, the longitudinal stress (a) diminishes 
to zero while shear stresses (x) rise to a maximum. As these ineffective 
lengths accumulate, loads continuously transfer from the weak (failed) 
fibers to the strong ones. If enough fiber failure sites are in close 
proximity to each other, the matrix begins to crack or de-bond from the 
fibers. When this occurs, the matrix can no longer transfer the load and 
the sample fails catastrophically. Thus composite strength can be 
modeled as a function of fiber strength, the ability of the matrix to transfer 
load, and the interface between the fiber and matrix. The mathematical 
model which predicts the probability of failure of the composite (in 
tension) from the statistics of fiber strength and the effective length has 
been established by Harlow and Phoenix [Ref 2 and 3]. 
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These three quanta are recognized as building blocks in composite 
strength determination. Matrix optimization is guided by this model. From 
a given assortment of matrix materials, one may determine which has the 
most suitable characteristics for stiffness and transferring shear stresses. 
Fiber strength statistics are determined through loading to destruction 
measurement. However, it is not feasible to directly measure and quantify 
the molecular level interface region. This can be accomplished by 
measuring the strength of a composite of known fiber and matrix, the 
fiber strength statistics and composite strength statistics can then be used 
simultaneously to solve for the ineffective length parameter. 



B. LIFE DETERMINATION OF FIBER REINFORCED 
COMPOSITES 

Decades of experience with homogeneous materials support the 
correlation between strength and life: that is, the larger the factor of 
safety for a given structure, the longer it can be expected to last. 

However, Coleman's theory of time dependence in the mechanical 
breakdown of fibers [Ref. 4] implies that the strength/life relation in load 
bearing fibers may not be so straightforward. Using the "weakest link 
hypothesis" first developed by Pierce in 1926, Coleman concluded that 
fiber life was instead a function of the statistical distribution of these 
"weak links" in time dependent strength. The time dependent strength of 
the fiber, termed breaking kinetics, is the extreme value of fiber life under 
a given load history. This is fundamental to the modeling of composite life. 
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Only when fiber life statistics are available can the model be applied to 
composite structures. Hence, the purpose for this work. 

The theory proposed is, the strength of a specific fiber can be related 
to the life of that same fiber. The dichotomy of course is, once a fiber has 
been broken in the strength test, another molecularly identical sample is 
not available for the life test. The alternative is to obtain a set of 
statistically identical fibers and test a portion of them in strength and 
another portion in life. 

This is being accomplished by using two different spools of AS-4 
graphite fiber, AS4-008 and AS4-019. A statistically meaningful number 
of samples have been taken from these two spools to estimate the 
respective parameters of the Weibull Distribution. The estimated 
parameters of these two different spools are different, suggesting the 
intrinsic strength of these two spools are from two different populations. 
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III. FILAMENT LIFE TESTING METHOD 



A. TEST SPECIMENS 

Test samples were single strands of Hercules Magnamite AS4 cut to 
two inch gauge lengths. Each sample consisted of one fiber filament 
mounted on a paper cardboard with a cutout for handling. The cardboard 
mount was severed with a hot wire prior to load application. Specific 
handling and mounting procedures are described in Appendix A and in 
Reference 5. Once mounted, the fiber diameters were measured using a 
laser diffraction device also following the procedures described in 
Reference 5. 

B. TEST EQUIPMENT DESCRIPTION 

The single filament samples were arranged in groups of 64 at four 
different sustained load levels by dead weights. The elapsed time 
between the startup of load and ending of stress-rupture at each station 
was triggered by an infra-red light switch. Test equipment is described in 
Appendix B. A polymethyl methacrylate (PMMA) cover was used to 
protect the samples from ambient air disturbances and ultraviolet (UV) 
radiation. It is known that UV light degrades the strength of some fibers; 
therefore, PMMA was chosen as a shield and is analyzed in Appendix C. 
Once loaded, the fibers were continuously monitored by a dedicated PC 
(personal) class computer. Software programs are reproduced in 
Appendix D. 
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C. TESTING PROCEDURE 

1. Strength Test 

Fiber strength data were acquired previously and are covered in 
References 5 through 7. 

2. Life Test 

During the diameter measurement process, every fifth fiber was 
tested to failure in strength using a calibrated load cell. This provided the 
link between the two groups needed to verify that the samples tested in 
life were statistically identical to those tested previously in strength. The 
actual loading of the samples on the board was carefully controlled and, 
as much as possible, were all done by the same operator. The loading 
procedure is also described in Appendix A. 

There were four levels of load set up on the board, 14.17 g. 
(Level 1), 13.15 g. (Level 2). and 12.20 g. (Level 3) and 10.96 g. (Level 
4). Each load level had 124 stations, 64 for each spool. This allowed a 
total of 512 life samples to be measured at sustained loads. Level 4 was 
planned but postponed. Instead those stations were also loaded to 14.17 
g. to gain additional measurements at that level. The sample board 
alternated rows with 008 and 019 fiber to give 4 rows of each. After 
each fiber was loaded, the date and time of the load were recorded in the 
lab notebook and also in a software spreadsheet. If a fiber failed on 
loading, its failure stress was recorded. After loading, each station was 
acdvely monitored by an AT computer which recorded stress failures as 
they occurred. Exact failure informadon was periodically retrieved and 
recorded in the notebook and spreadsheet. The spreadsheet 
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automatically computed realized data for stress ruptures, early censor 
(Type Iv) stress ruptures, fail on load stresses, and existing lifetimes for 
each fiber type at each load level. It would then arrange the desired data 
in a form easily transported to other programs for analysis or presentation. 
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IV. RESULTS 



A. TEST DATA 
1. Overview 

The data presented herein will be used to identify a suitable 
strength-life model for fiber filaments. A probabilistic formulation of the 
strength-life model is a joint distribution of the random variables, stress 
(applied), a, and elapsed time, t, denoted f(a,t). The parameters of this 
joint distribution are determined by three classes of realized data. The 
first is data realized in strength over negligible time. The second is data 
realized in life for a given stress history, in this experiment, constant, 
sustained stress. The third is data realized in strength given a survival in 
time history. The first two cases, that is, strength and life data, will be 
presented in this investigation. To examine variability in production, two 
spools of AS-4 graphite fibers were treated as two distinct populations, 
designated as spool 008 and 019. The division of fibers from these two 
spools is traced out in Figures 1 and 2. The entire population of either 
spool was denoted as {D}. This population is divided into strength test 
samples yielding data, {D s }, and life test samples yielding data, { } . The 

respective sub-division of these two data sets is described in the following 
two sections. 
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2. Strength Data 

The strength data set comes from work done previously using the 
same 008 and 019 fiber populations [Refs. 5 through 7], The sources of 
this data are summarized in Table 1. 

TABLE 1. SOURCES OF FIBER STRENGTH DATA 



DATA NAME REF. NO. 

SET 

{D . } Dr. Edward M. Wu 6 

Mr. Glen Nypiak 

{D s2 } Lt. David Bell, USN 5 

{D $ 3 } LCDR Carl Englebert, 7 

USN 

{D s 4 } Mr. Jim Nageotte 7 



3. Life Data 

Within the life data set, {D^ }, not all of the samples were realized 

in life. Some of the samples were tested to assure the life samples were 
statistically the same as the strength samples and are denoted as {D^ $ } 

(Figures 1 and 2). Some samples failed before they reached the target 
sustained load level. Thus, this was part of the life test subset but was 
realized in strength and is denoted as {D^-| s }, [D^h and {D^g}. Based 

on kinetic theory [Ref. 2], it is generally believed that the load supported 
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by a fiber increases as the loading rate increases; therefore, failures in 
strength would be hypothetically eliminated if the fibers could be loaded 
instantaneously. While this not physically possible, the data from samples 
which failed on loading can be accounted for by treating the life samples 
conditionally to strength greater than the sustained load. Data obtained 
from this experiment is summarized in Tables 2 through 4. Ruptures 
occurring on loading, prior to the desired load are identified as "Fail on 
Load". Stress ruptures occurring in time but caused by other than 
intrinsic stimulus internal to the sample are "Type Iv" censored. This data 
is statistically treated as greater than or equal to the realized life. The 
censorings are chance events as denoted by the variable time at censor. 
Stress ruptures in time are "Exact" and fibers still under load at the end of 
the experiment are "Type 1" censored. This data is statistically treated as 
greater than or equal to the time the data was interpreted. 



TABLE 2. RESULTS OF LOAD LEVEL 1 



Type of Failure 

Fail on Load {D-^ s } 



AS4-008 



AS4-019 



53 



57 



Type Iv Censor { D^ v } 



9 



7 



Exact {D{^ } 



15 



20 



Type I Censor { D^ c } 



116 



110 



TABLE 3. RESULTS OF LOAD LEVEL 2 



Type of Failure AS4-008 AS4-019 

Fail on Load {D t 2 S } 26 13 

Type Iv Censor 4 1 

{D t 2 v) 

Exact {D t2t } 9 6 

Type I Censor {D t2c } 62 63 



TABLE 4. RESULTS OF LOAD LEVEL 3 



Type of Failure AS4-008 AS4-019 

Fail on Load {D^ s } 10 11 

Type Iv Censor 1 3 

{D t 3v> 

Exact {D t3t } 3 1 

Type I Censor {D t3c } 62 63 



B. DATA ANALYSIS 

1. Analysis as a Function of Stress 
a. Overview 

Before proceeding with an analysis of life based on strength, 
it must first be confirmed that the two sets of data, {D s } and {D^ }, are 

indeed statistically from the same population. As can be seen in Figures 1 
and 2, each spool has yielded two sets of samples. The first set, {D s }, 
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was tested entirely in strength. Within the second set, {D^ }, every fifth 
fiber was tested in strength in order to link the two sample sets {D s } and 
{D t }. 

b. N on-Parametric Interpretation 

Non -parametric data analysis was conducted by plotting the 
cumulative failures in a weakest link (Weibull Probability) coordinate. 

This required the data to be ordered and a rank assigned. For this 
experiment, Expected Rank was used and is defined: the i th realization of 
the ordered xj is 

— for i = 1,2,3... N 
N + l 

and Xj is the ordered data (from weakest to strongest) in a set of N 
samples. An ambiguity arises when one considers the inclusion of 
censored data. Because the point is censored, its exact strength, and 
therefore its exact rank, cannot be precisely determined. Furthermore, 
because the data is ranked, this uncertainty perturbs the rank of each 
subsequent data point greater than the censored data. For this reason, 
censored data was omitted from the graphical representation. However, 
censored data did not affect the Likelihood calculations and was therefore 
included in the parametric interpretation. 

Once the data has been ranked, it is linearized by using the 
weakest link transformation. This method linearizes the ranked data and 
allows a visual observation of the lower tail characteristics. The weakest 
link transformation is defined as: 

F*(o) = ln(-ln(l-F(c))) 
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where F(o) is the empirical rank of the data in question. Recalling that 
each fiber is modeled as a chain of independent links, the strength of a 
given length of fiber must be the same as the strength of the weakest link 
in that chain, hence F*(g) can be further defined as the probability of 
failure of the weakest link of a given gauge length of fiber. 

Strength data, {D s }, from these two spools (008 and 019) is 

plotted in Figures 3 and 4. This set of fibers had previously been 
confirmed to have come from the same set of like materials and gauge 
lengths (Ref. 9). Figures 5 and 6 depict the F* plot of the set D ts overlaid 
with the set D s . The close proximity of the two sets of data indicates they 
too are of the same fiber and gauge length. A cursory run-test on the 
ranked set of merged data shown in Figures 7 and 8 shows acceptable 
levels of statistical clumping and completes the non-parametric 
identicalness test of the two sets. 

c. Parametric Interpretation 

A parametric analysis of the same data was also conducted 
by using the Weibull distribution for both strength and life: 



F(x) = 1 - exp 



Y 

x 



vP; 



where a and p are the shape and location parameters respectively. 

Given a data set, these parameters may be obtained by using the 
Maximum Likelihood Estimator (MLE) derived in Appendices E and F. 
The attribute of this method is there is no requirement to order the data 
which eliminates the process of ranking which allows the inclusion of any 
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censored data points. Also associated with the MLE is the concept of 
confidence interval. The confidence region of the parameters a and (3 
can be obtained from the likelihood contour plots. The confidence region 
is the confidence, stated as a percentage, that a selected a and (3 are the 
actual a and (3. Parameters chosen at a very low confidence interval 
have a low probability of matching the actual parameters, but when they 
do, they are highly accurate. A very high confidence interval, while 
allowing more room for error, is not very precise and is of little practical 
value. 

Figures 9 and 10 show a three dimensional representation of 
the Likelihood for the D s data set along with various confidence region 

contours. The 5 % and 1 % confidence region contours are shown in 
Figures 1 1 and 12. The shape parameter, a, and location parameter, (3, 
for this data set can be determined from Figures 3 and 4 graphically by 
the slope of the line of data points and the value of the X axis at F* = 0 
respectively. Alternately, the estimate can also be based on the calculus 
MLE derived in Appendix E. 

In order to show that {D s } is identical to {D-t s }, the likelihood 
ratio can be used. The likelihood ratio is graphically represented by the 
volume of the intersection of the two likelihoods at a given confidence 
interval, divided by the volume of the same confidence interval of the 
accepted baseline data set, { D s } or, symbolically 

V(D S )P1V(PJ 

V(D.) 
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: V(D S ) and V(D^ S ) are the volumes of the likelihood contours 

within the desired confidence interval. Figures 13 and 14 show overlay 
plots of the Likelihood contours of N s and N^ s at 5% and 1 % confidence 

intervals. By examining these figures, it can be determined that the 
likelihood ratios for 008 and 019 at the 5 % confidence interval are 
estimated to be .75 and .70 respectively, which is deemed to be 
acceptable for this confidence interval. 

The two sets of data can now be merged and the tested in 
life fibers are confirmed to be of the same population as the tested in 
strength fibers. As a last point of confirmation, Figures 15 and 16 overlay 
this merged data set {D s +D^ s } with the original set,{D $ } showing they 

are essentially identical, as do the contour overlay plots shown in Figures 
17 and 18. 

A distinction can now be made concerning those fibers 
which failed on load ({D^ ls }, { } - {Dt 3 s })• Although intended to be 

measured in time, these fibers actually failed in strength prior to reaching 
the target loads, , c^, and 03 respectively; therefore, this set of data can 

be considered as realized in strength to be included in the merged strength 
data set {D $ +D^ s }. The Expected Rank of the exact failed in strength 

data was presented as before with the exception that, the samples which 
reached the target load in life test were considered Type I strength 
censored at the stress rupture load. This is the equivalent to stating its 
strength would be greater than the stress rupture load. As with all non- 
parametric analyses, these censored points were not presented in the F* 
domain because of the ambiguity in rank. By treating censored data in 
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this manner, the graphical representation of the middle and upper tails 
remains unaffected, while the resolution of the lower tail is improved. Fail 
on loading data, {D^ s }, {D^g}, and {D^g}, is presented in Figures 19 
through 24 and can be analyzed in the same manner as the set {D^ s } with 
similar conclusions being drawn. Thus, the set of fibers tested in time can 
now be considered statistically identical to the set of fibers tested in 
stress. 

The final step in strength analysis was to analyze the entire 
population of strength data and determine a shape (a) and location ((3) 
parameter for the Weibull Probability of Failure Function which would be 
necessary to determine the Likelihood in Life. Results of non-parametric 
analysis of the set {D} in strength are shown in Figures 25 and 26. 
Because of the data censored above 12.20, 13.15, and 14.17 grams, a 
value for a and (3 can only be estimated from the set {D s + D^ s }; 

therefore, a parametric analysis in the form of the Likelihood contours of 
the set {D} is shown in Figures 27 and 28. 

2. Analysis as a Function of Time 

Since the life of a fiber is dependent on the fiber's surviving the 
loading process, ruptures during loading must be considered in 
determining the Likelihood of the fiber in life. The reliability of the fiber in 
strength will be necessary to determine its reliability in life. Therefore, 
equation (1) in Appendix E must be modified to include the fibers 
reliability in strength: 
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R(o) = exp' 
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and the Likelihood in the time domain becomes: 
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where o c is the target stress and k is the number of left censored fibers. 

By determining estimates of ag and (3a from Figure 27 for Spool 008 and 
evaluating the Likelihood contour over a set range of oq and (3^, the 

Likelihoods for each of the three load levels can be plotted and are shown 
in Figures 29 through 31 which were developed using an estimate of 5.0 
for a a and 18.0 for (3a. A similar procedure was followed for Spool 019. 

Although not investigated during this analysis, an alternative 
method of calculating the Likelihood would be to evaluate the fail on 
loading data as censored in time rather than stress. In order to do this, 
one must first assign a time, tg, at which the sample is considered to have 

attained the target stress. All samples which fail on loading can be 
considered to have a life less than or equal to tg, or: 
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F(t) = 



1 - expj 



a , 

ft ^ 

o 



P 



V^t J 



Under this interpretation, the Likelihood can be calculated using: 



F(t) = 



1 - exp 



a i 

ft ^ 



P, 



v^t j 



-i k 



„ m a m i , r —a i 

a t 1 P t nx i 1 ex p{- P, x “'} 

i = l 



fl exp 
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3. Construction of the Strength-Life Model 

For a life distribution given a sustained stress level, Gj, the Weibull 

model used was: 



F(tla = a.) = 1 - exp' 



f _L^ 

P, 



Using this model, the Likelihood contours for the life data realized up to 
the current time were plotted and are shown in Figures 29 through 31.The 
flatness of the contours is the result of an extremely broad range of |3 
which is caused by the elapsed time (at the time of data interpretation) 
being much smaller than the underlying location parameter of life, (3. 
Therefore, only very tentative estimates of (3 are possible at this time. 

The shape and location parameters are visually estimated to be: 
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Stress Level, a. 


6t, 


P, 


14.17 g. 


.15 


l.OE+11 


13.15 g. 


.15 


1.0E+12 


12.20 g. 


.15 


5.0E+13 



The estimated life at different levels of sustained stress is shown 
in Figure 32 evaluated at three levels of probability of failure, 0.001, 0.01, 
and 0.1. The figure can be considered as a graphical representation of the 
Strength-Life Model for AS4 graphite under sustained stress. The 
expected life of the fiber, evaluated at the desired probability of failure at a 
given stress level, can be obtained directly from the model. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The Strength-Life Model, once verified, can be a useful and important 
tool for predicung fiber reliability using strength statistics. In design, it 
allows accurate determination of maximum design stress for an expected 
life. In maintenance and repair, it facilitates the calculation of the 
remaining life of an existing structure. It also predicts stress levels at 
which an existing structure may be operated that can increase or 
decrease its remaining life. 

The following items are recommended for follow-on research:. 

1. Application and verification of the model using larger scale 
composite samples. 

2. Mathematical formulation to predict life based solely on strength. 

3. The life test be continued to improve the estimation of the 
parameters of the model. 
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APPENDIX A 



PROCEDURES. 

Fiber Handling. Fiber handling procedures are specific and are in 
effect continuously. 

1. Unless required for the immediate task at hand, all fibers and 
filaments shall be stored in dark, secure containers. 

2. Fibers and filaments are fragile. Any fiber or filament that is 
dropped or otherwise subjected to extreme stress, shall be discarded. 

3. Whenever possible, all samples should be handled by one end 
only so as not to exert any tension on the sample. 

4. Use of fluorescent lights shall be kept to a minimum. 

Filament Mounting. Filaments were prepared according to the 
following guidelines: 

1. Strip off a random length of fiber from the spool, sever, and 
discard. 

2. Cut the next 12-15 cm. of fiber from the spool and tape one end 
to the glass plate. 

3. Pour a small puddle of ethyl alcohol on the glass plate to float the 
free end of the bundle. The filaments will begin to float apart. Gently 
fan out the filaments floating in the alcohol and let stand until all the 
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alcohol has evaporated. 

4. Write the correct serial numbers on the carriers to be loaded. 

5. Using a small piece of transparent tape on the end of a modelling 
knife, lift the free end of one filament and slide an empty carrier under 
it 

6. While maintaining filament alignment with the two small holes in 
the carrier, lower the raised filament end and tape it to the glass 
plate. 

7. Use two small pieces of tape to secure the filament to the 
carrier. 

8. After the filament is secured between the two pieces of tape, 
carefully cut the filament to free it and the carrier from the glass plate 
and place the filament/carrier assembly on the cold curing plate. 

9. Repeat steps 5-8 eight more times. After the ninth filament has 
been mounted, discard the remains of the bundle. 

NOTE 

In order to maintain standard curing times and temperatures, do 
not activate the heat element in the hot plate until the plate is full and 
do not remove the samples until the plate has been turned off and has 
cooled. 

10. Activate the hot plate with a setting of 200° F. 

1 1. Thoroughly mix a small amount of epoxy (Devcon 2 Ton). 

12. Reduce hot plate temperature to 150° F. 
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13. Place one dot of epoxy at each end of the elongated hole in 
each carrier to permanently secure both ends of the filaments to their 
respective carriers. Also, place a dot on the outboard ends of the two 
small holes to re-enforce them. 

14. When the last filament has been epoxied, turn off the hot plate 
and let cool. 

15. Store mounted filaments in designated containers. 

Filament Loading. 

Loader Start up 

1. Turn on power supplies for load cell and elevator. 

2. Turn on HP-85 and attached Data Acquisition Unit. 

3. Insert program tape into HP-85 and type LOAD "QUIK4". 

4. Type RUN. 

5. Select "k2" (INPUT) key and enter the date. 

6. With load cell assembled and unladen, select "kl" (ADJ B) key. 
This sets the bias to 0. To check, select "k3" (WEIGH). If a weight 
within .003 g of O.Og is not obtained, re-select ADJ B until satisfactory. 
Accurate readings should also be verified using calibrated weights. 

Filament Loading 

1. Remove Plexiglass cover from sample rack. 
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2. Without disturbing adjacent filaments, lower stabilizer bolts as far as 
possible on stations to be loaded. 

3. Adjust the spring so that the flag will operate the optical trigger 
without rubbing against it. 

4. Remove the weight from the station to be loaded. 

5. Using tweezers, grasp the desired filament by one end of it’s carrier 
and hang it on the hook attached to the flag. 

6. Using tweezers, hang the weight on the bottom hole of the carrier. 

7. Verify the flag does not rub against the trigger. 

8. Lower the load cell elevator far enough that it will fit onto the rack 
without disturbing the hung filament. 

9. Without anything touching the top of the load cell, press "ADJ B" on 
the HP-85. 

10. When the screen has stopped scrolling, raise the elevator until the 
compliance spring is completely retracted and the paper carrier is 
carrying exactly zero load. 

1 1. Verify that the PC records the station as unloaded as the flag 
activates the optical trigger. 

12. Select "LOAD" on the HP-85 and answer the cues: 

"ENTER THE SAMPLE #" 

"ENTER THE STATION #" 

13. Select "CONT" 

14. Verify the weight of the vial shown by the computer is correct for 
the station being loaded. 

15. Ensure the direction switch on the elevator is in the "down" 
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position and the power switch "off'. 

16. Check the setting on the hot wire power supply and turn it on. 

The wire should glow a dull red, not bright. 

17. Without touching the filament or disturbing adjacent filaments, 
carefully burn through both sides of the paper carrier at the bottom of 
the slot. 

18. Turn off the hot wire and return it to it's holder. 

19. Select "CONT" on the HP-85 and move the elevator power 
switch to "on". Monitor the elevator as it lowers to ensure it does not 
bottom out prior to being unloaded. 

20. Verify the PC records the station being loaded as the flag blocks 
the optical trigger. 

NOTE. If the filament breaks before it is completely loaded, stop 
the elevator, record the mass displayed by the HP-85, and start again 
from step 24. 

21. When the weight hangs free, the HP-85 will beep twice. When 
this happens, stop the elevator. After two minutes, the HP-85 will beep 
twice again, allowing you to continue. 

22. Lower the elevator until it can be withdrawn without disturbing the 
filament. 

23. Once the elevator has been withdrawn, carefully, without 
disturbing the filament, adjust the stabilizer screw so that it is 
centered on the vial without touching it and has a vial/stabilizer gap of 
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not more than 1 mm. 

24. Load a sheet of paper into the plotter and select "shift k2" (PLOT) 
on the HP-85, answering the cues as required. 

25. Record filament load information in Fiber Life Data Book and file 
load plot in the plot binder. 

26. Repeat from step 2 as necessary. 

27. When all desired stations have been loaded, verify all stabilizer 
bolts are positioned correctly and replace Plexiglass cover. 

28. Verify PC is reporting the correct status for the stations loaded. 
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APPENDIX B 



EQUIPMENT DESCRIPTION 

Data for the experiment was collected by suspending fiber samples on 
the individual stations (512 total) of a large sample board. Load and fail 
status was actively monitored by individual infrared switchesthrough 
ribbon wire which fed the data to an interface box communicating with a 
PC- AT. (Figure B-l). 

The sample board was made up of 512 stations arranged in 8 rows 
and 64 columns. Each station consisted of a mounting point, a high 
compliance spring (k=0.635 mm/g), an optical switch and trigger, a 
filament sample, a weight, and an adjustable stabilizer screw. The board 
was constructed of aluminum channel stock and fiberboard. When not 
being loaded, the stations were covered with a protective sheet of 1/4" 
PMMA. 

The trigger was a wire paper clip with the inner loop removed. The 
top half of this outer loop was then wrapped in black vinyl electricians 
tape. The fiber sample, mounted in its paper carrier, was hung on this 
trigger which in turn hung from the, bottom of the spring. With the weight 
attached to the bottom of the fiber carrier, the black tape blocked the 
optical switch. When the filament broke, the weight dropped away while 
the spring retracted the trigger allowing light to pass through the paper 
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clip thus closing the optical switch, which was detected and recorded by 
the PC. 

The weights were plastic vials filled with lead shot and were loaded to 
3 weight levels: 12.20, 13.15, and 14.17 g. Two rows of AS-4 graphite, 
batch no. 008, and two rows of AS-4 batch no. 019 were loaded at 14.17 
g., one row of 008 and one row of 019 were loaded at 13.15, and one row 
of 008 and one row of 019 were loaded at 12.20 g. It was decided to 
double the number of high weighted samples in order to obtain more data 
points prior to censoring the experiment. 

Filaments were permanently mounted in paper carriers to facilitate a 
reliable and convenient method of storage and handling. Once assembled 
in it’s proper station, the paper carrier was severed using a hot wire 
instead of scissors to reduce the chance of sending a shock through the 
filament. [Ref. 3] 

Each of the 512 optical switches was connected to a custom designed 
interface box. This box contained a power supply for the optical switches 
as well as a multiplexing circuit which fed status information to a 
MetraByte I/O card in the PC. The board was divided into 64 8-bit 
words and software was written to cycle through each of the words and 
to report any changes in the value of the word. For example, if stations 1- 
8 are all loaded the PC reads a value of 0000 0000. If station 2 fails, the 
next time this byte is read, it would show a value of 0000 0010. The PC 
would then record the failure and continue monitoring the board. A back- 
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up power supply was connected to the interface box and PC in case of a 
filament failure during a power outage. PC software is reproduced in 
Appendix D. 

Start times for each filament were also recorded by the PC when a bit 
change from 1 to 0 was detected. However, consistently accurate start 
times were highly dependent on a function of individual spring compliance, 
precise trigger construction , loader motor speed, and masses of the 
weights. In order to eliminate this variability, the loader software used to 
program the HP-85 also incorporated a timing mechanism. As the load 
cell was lowered, the HP-85 monitored the weight on the cell. When zero 
weight was detected, a beep sounded and a two minute timer was 
started. The load cell was left under the sample until expiration of the two 
minute time period. Should a sample have failed during the first two 
minutes of life, its lifeume would be recorded to the nearest 1/10 second. 
After two minutes, such precise record keeping was not necessary and 
timing was turned over to the PC. 



The HP-85 was programmed to read the load cell approximately 5 
times per second during the loading sequence and could therefore provide 
accurate data on the loads as they were applied. After loading, this 
information was sent to a plotter and a record of each load made for 
analysis. If a filament ruptured any' time during the load sequence, up to 
the expiration of the two minute time period, the HP-85 would record the 
load on the filament at the time of failure or, if already loaded, the life of 
the filament up to the first two minutes. 
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Filaments were loaded by using a 150 g. capacity Sensotec load cell 
mounted to an electrically controlled, hydraulically actuated elevator. The 
electric motor, Western Gear model P5B24R3 incorporated a built in 
reduction gear and was run at 7.5 volts giving a start to finish load run of 
about 30 seconds yielding 170-220 data points as the load on the filament 
increased. The motor was mechanically connected to a disposable 5 cc 
hypodermic syringe which acted as the master cylinder. The slave 
cylinder was an identical syringe with the load cell mounted on a platform 
secured to the plunger. This arrangement was in turn mounted to a 
Plexiglass platform which provided a stable base during loading. The 
slave and master were connected with 3 feet of 1/8" ID vinyl tubing; the 
hydraulic fluid was tap water. 
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Figure B-l. Diagram of Fiber Life Testing Apparatus. 
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APPENDIX C 



ULTRAVIOLET (UV) LIGHT SHIELDING 

The life of certain fibers, (Kevlar for example) are known to be 
affected by exposure to UV radiation. To assure these life test 
experiments were performed in an inert environment, an investigation into 
the presence of UV light was conducted. 

The sample rack was assembled in a secure laboratory environment 
on a basement level with windows facing Northeast. The rack itself was 
mounted on the Southeast wall and therefore, could not receive any direct 
sunlight. Windows were covered by opaque commercial grade drapery. 
Room lighting was fluorescent overhead with plastic covers. The sample 
racks were covered with 1/4 inch Polymethyl methacrylate (PMMA) 
Plexiglas. Light from the two sources was filtered twice before reaching 
the samples. The indirect sunlight passed through standard grade soda 
lime window glass and then through the PMMA rack covers. Overhead 
lighting was filtered through a plastic optical diffuser and then through the 
PMMA rack cover. 

Tests for the intensity of 300-400 nm UV light w-ere conducted using a 
100 Watt mercury lamp at a distance of 30 ft. from the spectrometer. All 
measurements were made at the surface of the sample material. Results 
of the spectrum analysis are shown in Figures C-l through C-3. To 
explore materials for improved UV shielding, analysis was also done on 
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Polycarbonate (Lexan) and a UV filtering grade of PMMA (Acrylite OP- 
2). Results from these two materials are shown in Figures C-4 and C-5. 
The improvement in shielding appeared to be a small range of an 
additional 25 nm to that already provided by the PMMA and was 
therefore not deemed to be significant. 

Based on these measurements, the amount of UV light reaching the 
actual fiber samples was judged to be insufficient to affect life testing. 
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Transmittance: Soda Lime Glass 
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Figure C-l. Results of UV Transmittance test for Soda Lime Glass. 
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WAVELENGTH (nm) 



Transmittance: PMMA (Plexiglass) 
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Figure C-2. Results of UV Transmittance test for PMMA. 



36 



WAVELENGTH (nm) 



Transmittance: Flourescent Lighting Cover 
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Figure C-3. Results of UV Transmittance test for Installed Flourescent 
Lighting Cover. 
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WAVELENGTH (nm) 



Transmittance: Polycarbonate (Lexan) 




Figure C-4. Results of UV Transmittance test for Polycarbonate. 
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| Transmittance: UV Filtering PMMA (Acrylite OP-2) 
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Figure C-5. Results of UV Transmittance test for UV Filtering PMMA. 
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APPENDIX D 



MAIN PROGRAM: 

' ** Filamod7 ** 

* Jim Nageotte Nov 13, 1990 



Start: 

DIM newrd% (64) , 01dRd%(64), row%(8), 
Status$ (2) , StarTime$(512) 


* Dimension Variables 

* 



row%(l) = 1: row%(2) = 2: row%(3) =4: * Declare constants 

row% ( 4 ) = 8: row% (5) = 16: row%(6) = 32: 
row%(7) = 64: row%(8) = 128 
a% = 0 



OUT &H303, &H10 * 


* Reset I/O Card 


TYPE S tatus2 

OldRd AS INTEGER 
END TYPE 

DIM Stat AS Status2 


* Declare file type 

* for status file 


TYPE Times 

StarTime AS STRING * 16 
END TYPE 

DIM Tim AS Times 


* Declare file type 

* for Start times file 


CLS 

PRINT "Ver 11.13.90" 

GOSUB Gettime 

PRINT "Power up at TIME$ ; " on " ; 

DATES 

PrintStat = INP(&H3BD) • 

IF PrintStat = 223 THEN * 

LPRINT "Power up at " ? TIMES? 


* Print start up message 

* & current date on screen. 

* Read printer status. 

* 223 indicates that the 

* printer is ready. 


" on " ? DATES 
LPRINT : LPRINT 
END IF 


* Print start up message 

* only if printer is ready. 


OPEN "Status. Dat" FOR RANDOM 
AS #1 LEN = 2 * 

FOR idx = 1 TO 64 1 
GET #1, idx, Stat 
01dRd% ( idx) = Stat. OldRd 
NEXT idx 


* Open Status File for 

* random access. 

* Read previous status 

* from disk file. 


CLOSE #1 


* Close Status File. 


OPEN "Times. Dat" FOR RANDOM AS #2 
LEN = 16 ' 

FOR i = 1 TO 512 
GET # 2, i, Tim 
StarTime$(i) = Tim. StarTime 
NEXT i 


* Open Start Times File 

* and read start times 

* for all 512 stations. 
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CLOSE #2 



* Close Start times file. 



GOSUB Checksw 
FOR i = 1 TO 64 



* Read 8 bit status words. 

* until something changes. 



IF newrd% ( i) = OldRd%(i) THEN GOTO Same 
a0% = 01dRd% ( i) : al% = newrd% (i) 

FOR j = 1 TO 8: jdx% = 9 - j 
b0% = FIX (a0% / row% ( jdx%) ) 
bl% = FIX (al% / row% ( jdx% ) ) 

IF b0% = bl% THEN GOTO Nextbit 
Station = (i - 1) * 8 + jdx% 

IF bl% < 0 OR bl% > 1 THEN 

PRINT "Error in Status bl%" 



* Compares Status read from 

* disk to that read from 

* the switches. If they are 

* not equal it looks more 

* closely at the switches. 

* 0 or 1 are the only valid 

* values for bl%. 



BEEP 
END IF 

IF bl% = 1 THEN 
PRINT "Station #" ; Station; 
"Loaded RdTime$ ; " on " ; 

RdDate$ 

GOSUB Store 

PrintStat = INP(&H3BD) 

IF PrintStat = 223 THEN 
LPRINT "Station #"? 
Station; " Loaded "; 
LPRINT RdTime$; " on "; 
RdDate$ 

END IF 
BEEP 
END IF 



*"1" indicates "Loaded" 

* Print loaded message to 

* screen. 

* Store new Start time to 

* disk file. 

* Check printer Status. 

* If printer is ready 

* then print loaded 

* message to printer. 



IF bl% = 0 THEN 
GOSUB Elapsed 

PRINT; PRINT "Station Station; 

" "; StarTime$ (Station) ; "to" ; PRINT 

LEFT$ (RdTime$ , 5); " " ;RdDate$ ; PRINT 
" " ; YEARS; "Years, "; DayT; "Days " 
PRINT Min; " Minutes": PRINT 
PrintStat = INP(&H3BD) 

IF PrintStat = 223 THEN 

LPRINT " "; LPRINT "Station 
Station;" "; 

StarTime$ (Station) ; " to "; 

LEFT$ (RdTime$ , 5);" RdDate$ 
LPRINT " "; YEARS; "Years, "; DayT; 
LPRINT Min; "Minutes". 

LPRINT " » 

END IF 



* "0" indicates failed. 

* Find elapsed time. 

* Print failure information 

* to screen. 



* Check printer Status. 

* If printer is ready, 

* then print failure 

* information to printer. 



"Days 



T < - 
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END IF 



Nextbit : 

aO% = aO% - bO% * row%(jdx%): al% = al% - bl% * row%(jdx%) 

NEXT* j 

01dRd% ( i) = newrd%(i) 

OPEN "Status . dat" FOR RANDOM AS 

#1 LEN = 2 * Update disk Status file. 

Stat.OldRd = OldRd% ( i) 

PUT #1, i, Stat 
CLOSE 

Same : 

NEXT i 
LOOP 

END * End of loop that compares 

* old status to new and end 

* of main program. 

SUBROUTINES : 



Check Switch compares the current status to the status read off the disk 
file. It loops until there is a change. It then sets a flag and 
returns with the current time and date to start checking each of the 512 
stations and determine what type of change has taken place. By checking 
the status in 8 bit words and in such a small loop, we keep the time 
required to check all switches and the response time to a minimum. 



Checksw: 

DO 

GOSUB Gettime 

a% = 0: idx = 0: RdTime$ = TIMES 
RdDateS = DATES 

FOR i = 1 TO 8: OUT &H301, row%(i) 
FOR j = 0 TO 7: OUT &H302, j 
idx = idx + 1 
newrd% ( idx) = INP(&H300) 

IF newrd%(idx) <> 01dRd%(idx) 
THEN a% = 1 

NEXT j 
NEXT i 

LOOP UNTIL a% - 1 ‘ 

BEEP 

RETURN 



* Get current time. 

* Get current date. 

* Select Row to read. 

* Select word to read. 



* Set Flag if Status has 

* changed. 



* Loop until Status has 

* changed. 



42 



Elapsed Time calculates the time between the Start time read from disk 
and the current time. It calls two other sub routines to. do this, Julian 
and Minutes . Julian is called twice, the first time it converts the 
Start time to a Julian date and the second time it converts the current 
date to a Julian date. The conversion is necessary to allow subtraction 
to be performed on the dates. Having arrived at the number of Years and 
Days before the fiber ruptured, we call the subroutine Minutes to 
subtract the times involved and return with the correct number of 
minutes. An adjustment to the number of days may be made based on the 
number of minutes. 



Elapsed : 

ConvDate$=RIGHT$ (StarTime$ (Station) , 10) 

GOSUB Julian 

StartYr * Year 

StartDy = Days 

ConvDate$ = RdDate$ 



GOSUB Julian 

StopYr = Year 

StopDy = Days 

YEARS = StopYr - StartYr 

DayT = StopDy - StartDy 

IF DayT < 0 THEN 

DayT = DayT + 365 
YEARS = YEARS - 1 
END IF 

GOSUB Minutes 
IF Min < 0 THEN 

Min = Min + 1440 
DayT = DayT - 1 
END IF 
RETURN 



* Select Start Date. 

* Convert to julian date. 



* Select Current Date. 

* Convert to julian date. 



* Subtract years. 

* Subtract days. 



* Calculate minutes. 



Julian: 

Mon = VAL(LEFT$ (ConvDate$ , 2)) * Get month, day and year 

Day = VAL (MID$ ( ConvDate$ , 4,2)) * from the date string. 

Year = VAL(RIGHT$ (ConvDate$ , 2)) 

SELECT CASE Mon 
CASE 1 

Days = Day 
CASE 2 

Days = Day +31 
CASE 3 

Days = Day + 59 
CASE 4 

Days * Day +90 * Calculate total days. 

CASE 5 

Days = Day 4- 120 
CASE 6 

Days = Day + 151 
CASE 7 

Days = Day + 181 
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CASE 


8 








CASE 


Days 

9 


= Day 


+ 


212 


CASE 


Days 

10 


= Day 


+ 


243 


CASE 


Days 

11 


= Day 


+ 


273 


CASE 


Days 

ELSE 


ll 

D 


+ 


304 




Days 


= Day 


+ 


334 



END SELECT 
RETURN 

Store simply opens the Start Time file on disk, a random access file, 
and stores only the new time without reading or re-writing the 
unaffected times. The file is then closed and the routine exited. 

Store: 

OPEN "Times. Dat" FOR RANDOM AS #2 LEN = 16 
Tim. StarTime = LEFT$ (RdTime$ , 5) + " " + RdDate$ 

PUT #2, Station, Tim 
CLOSE H 2 

StarTime$ (Station) =■ Tim. StarTime 
RETURN 



Minutes recovers the values of Minutes and Hours from the time strings 
and converts them to Minutes only for ease of mathematical operations 
in the time keeping. 

Minutes : 

HrSt$ = LEFT$ (StarTime$ (Station) , 2) 

Min$ = MID$ (StarTime$ (Station) , 4, 2) 

Hrs = VAL(HrSt$) 

MinSt = VAL(Min$) 

MinS = Hrs * 60 + MinSt 
Hrf $ = LEFT$(RdTime$, 2) 

MinF$ = MID$(RdTime$, 4, 2) 

Hrf = VAL(Hrf $) 

MinFn = VAL(MinF$) 

MinF = Hrf * 60 + MinFn 
Min = MinF - MinS 
RETURN 



Gettime reads I/O ports to read the time directly from the system clock. 
It then converts these values to a string that resembles the "TIME$" 
function. The program then uses this string called ,f Rtime$" for Real 
Time String. 
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Gettime: 

Sec$ = LTRIM (STR$ ( INP ( &H2C0) AND &HF) ) 
Secs$ = LTRIM (STR$ (INP ( &H2C1 ) AND &HF) ) 
Minute $ = LTRIM ( STR$ (INP (&H2C2) AND &HF) ) 
Minutes$ = LTRIM ( STR$ (INP (&H2C3) AND &HF) ) 
Hour$ = LTRIM (STR$( INP (&H2C4) AND &HF) ) 
Hours $ = LTRIM (STR$( INP (&H2C5) AND &HF) ) 
Rtime$ = "00:00:00" 

MID$(Rtime$, 1) = Hours$ 

MID$ (Rtime$ , 2) = Hour$ 

MID$ (Rtime$ , 4) = Minutes$ 

MID$ (Rtime$, 5) = Minute$ 

MID$(Rtime$, 7) = Secs$ 

MID$ (Rtiroe$, 8) = Sec$ 

RETURN 
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Figure D-l. Flowchart for Life Test Monitoring Program. 
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APPENDIX E 



I. LIKELIHOOD DERIVATION FOR EXACT DATA 



For Exact Data: 



L exact = f ( x i’ a ’ P) = d/dx F(Xj, a, (3) 



:F(x,a,p) =1 




Taking the derivative with respect to x: 



f (x) = |E = -eK| 



al 



( 1 a-W 

“a-CXX 

U J 



a/_x 

pU 
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The Likelihood for all data, x, , = 



-( x i) = .n f (\- ^ P) 
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Rearranging, 
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Let A = InL = In f(Xj,a,P) 



= ln^a"f 



ln|a n p' nu flx i (a - ,) eKp{-p' a ix“}j 



= n lna-nalnp + Xln 



i=i 



ex P {- rix“}] 



= n Ina - not ln|3 + (a - 1)1 n X x t + Xlnfeup-^- )' Ix“ 




= n lna-nalnp + (cx- 



n 

v 



Din 2. x 

i=1 1 



-a n 
-p I 

i=1 



x a 

i 



(3) 



To find the maximum Likelihood, (L max ), set the derivative of A with 
respect to a and (3. to zero: 



Taking the partial derivative of eq. (3) with respect to a and (3 and 
setting it to zero yields: 



n In p+ Z^l nx . + (3 lnpX^x“-p 



-a 



a 



n 

X x: 

i=l 1 



Inx. =0 

i 



(4) 



0 A _ _ not Ct y a _ n 

dp P ' 



(5) 



Rearranging eq. (5): 



n a a 



p 



lx 



. a 



1 n 

n =Ay x a 
p i=l 
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1=1 J 



( 6 ) 



(7) 
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Now replace the n in the 2nd term of eq.(4), with eq. (6): 



n_ 

a 



f A ^ ^ —Cl ^ nr ^ 

-4r Ix a lnp + Sin x + p lnpSx a -p Ix c 
^ P i= 1 1 J i= 1 1 i =1 1 i =1 1 



and cancel with the 4th term: 

77 + ilnx. -p" a ix a lnx = 0 

a >=i 1 i= i 1 

From eq.(6), substitute 

a n 

np = Ix a 
1= i 1 

and divide both sides by n yielding: 



1 

a 



1 n 

+ F lll)x i 
i =1 



I x a Inx 
1 

n 

x x a 

, i 



Inx. = 0 

1 
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Rearranging to set the right side equal to a: 



1 

a 



£x a lnx 
* 1 1 



Lx a 
1=1 1 



i n 

-Ti-Zlnx. 

i =1 



,or 



a = 



Sx a lnx 

i=i 1 



l i n 

s— -F.V’ 1 . 

Ix u 
1=1 1 



( 8 ) 
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APPENDIX F 



I. LIKELIHOOD DERIVATION FOR EXACT DATA WITH 
TYPE 1 AND TYPE I-V CENSORING 

For Exact data which has both Type I and Type I-v (Variable: Cause of 
failure not quantifiable in this experiment) censoring: 

Likelihood. = []~[f(x i ,a,p)]]^[[l - F(x,,a,p)][l - (n ~ m) 



Exact Type I-v 



Type I 






n =The total numberof fibers 

A 

m =The total numberof fibers that failed 
m! = The numberof fibers failing in life 



m - m ( =The numberof fibers subject to Type I - v 



censoring 



Therefore, 



L = 



ni| 

a m -p--rixr i exp{-rxr} 

1 =) 



11 . 

n-P 



i=m, 






H J 



exps- 






-i(n-m) 



\P J 



( 1 ) 



Let A s In L 



= In 



a 



-m | a m i 



p 1 n x“ 1 eKp{- p~V} 



In 



m 

n e«p 



i =m 



i 



f £ oT 




1 (A] 


“11 


(n-m)" 


pj I 


+ In 


i eKp iiTi 







m m 

a 



= m Ina - m a Inp + (a - 1 ) Z Inx -p Xx a - 

i=i 1 1=1 1 

a 

/ V 

m 



P Xx -(n-m) 



x . 
j 

P J 



( 2 ) 
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To find the maximum likelihood of L, 



Set #■=# = * 

da ap 

Taking the partial derivative of eq.(2) with respect to a yields: 



m , 



|^ = -^ L -m 1 lnp+ I Inx -] p a Ix“lnx ; - p"“lnpix“ 



i= 1 



m , 



i= 1 



i=1 



m 



P" a I x^lnx . - p "In p I £ 



m 



i =m 



i= m 



yv a 
i 



- (n - m) 



fx.\ fx.\ 

i) “(t) 



(3) 



Taking the partial derivative of eq.(2) with respect to (3 yields: 



m , 



^ =- m.ap" 1 + ap" a_1 x'x^ap -0 ' 1 I £“ + ap"“ -, (n - m)x a = 0 
op 1 ,-1 1 ,_ m 1 J 



1=1 



Multiplying through by 



. a+1 



a 



m j 

-m,p a +Ix a + I x “ + (n - m)x a = 0 
1 .=i 1 um, 1 J 
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Now divide by mj and rearrange to isolate (3: 






m 



Ix“ + 
i=i 1 



m 

I 



x“+ (n - m)x“ 




To solve for a, first multiply eq. (3) by (3 a and rearrange: 



a 



- m 





Now gather the ln(3 terms. 



m 



rn , 



/ 



m , 



p — — mlnp+XInx + In P| X x u + Xx +(n-m)x a 
V i=J ') [i=l 1 i= m| 1 \ 

m i m 

X x a lnx + X xlnx +(n - m)x a ln 

.= i 1 ' 



1 1 

1=m , 



m 
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Divide both sides by mj and rearrange to bring a to the left hand side: 



1 

a 




Inp + pj^-XInx 



i i=i 



(5) 



Obtain from eq. (4), 

m , 

f 1 m 

m R = Xx H + X x + (n -m)x u (6) 

1 ,=1 1 i= m( 1 J 



and substitute for (a) in eq. (5): 



1 1 

a m 
m R 



rn , 



m 



Xx“lnx.+ X x Inx +(n-m)x (I lnx 



. , 1 1 1 1 
1=1 i =m , 



m 



[m.p-Jnp 



!• + 



i i 

InB + — X Inx 
m n=i 1 
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Cancel the lnp terms and invert both sides: 




Finally, eliminate P by substituting eq. (6) for the denominator: 
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APPENDIX G 




Figure 1 . Relationship of Fiber sample sets from Spool AS4-008 
population. 




Figure 2. Relationship of Fiber sample sets from Spool AS4-019 
population. 
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Figure 3. Expected Rank of 008 strength set {D s } 
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Figure 4. Expected Rank of 019 strength set {D s } 
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Figure 5. Expected Rank of 008 strength sets {D s } and {D ts } 
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Figure 6. Expected Rank of 019 strength sets {D s } and {D ts } 
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Strength - s [gm] 





A 


B 


C 


1 


AS4-008 






2 








3 








4 








5 








6 


xi (Ds|Ns) 


xi (Dts|Nts) 


F*({Ds+Dts}|N 


7 




5.909 


-6.277309565 


8 


6.393 




-5.382913936 


9 


6.666 




-4.917448539 


10 




6.711 


-4.600324093 


11 


7.12 




-4.359337455 


12 


7.997 




-4.164810785 


13 


8.017 




-4.001618953 


14 




8.54 


-3.860998686 


15 


8.9 




-3.737415596 


16 




9.085 


-3.627150043 


17 




9.364 


-3.52758132 


18 


9.58 




-3.436793307 


19 


10.037 


— 


-3.35334254 


20 


10.04 


[ 


-3.276114499 


21 




10.047 


-3.204230672 


22 




1 0.1 36 1 


-3.136986293 


23 


1 0.234 




-3.073807343 


24 


1 0.263 




-3.014220084 


25 




10.551 1 


-2.957828979 


26 




10.564 


-2.904300374 


27 


10.588 




-2.853350225 


28 


1 0.592 




-2.804734725 


29 


10.654 




-2.758243047 


30 


10.776 




-2.713691648 


31 




10.909 


-2.670919757 


32 




10.929 


-2.629785756 


33 


10.973 




-2.590164254 


34 




10.987 


-2.551943695 


35 


10.998 




-2.515024401 


36 


11.169 




-2.479316935 


37 


11.199 




-2.444740751 


38 


1 1 .279 




-2.41122305 



Figure 7a. Expected Rank data of merged 008 strength set {D s +D ts } 
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39 


11.378 




- 2.378697815 
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11.4 




- 2.347104998 
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11.417 




- 2.316389819 
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- 2.286502166 
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11.483 




- 2.25739608 


44 


11.485 
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45 




11.548 


- 2.20136292 


46 




1 1 .598 
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47 
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1 1 .808 
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1 1 .869 
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51 


12.043 
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53 


12.141 




- 2.00143457 


54 
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55 


12.18 
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56 
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57 
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59 
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63 
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64 
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65 
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- 1.56501914 



Figure 7b. Expected Rank data of merged 008 strength set {D S +D( S } 

(cont.) 
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Figure 7c. Expected Rank data of merged 008 strength set {D s +D ts } 

(corn.) 
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Figure 7d. Expected Rank data of merged 008 strength set {D s +D ts } 

(corn.) 
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Figure 7e. Expected Rank data of merged 008 strength set {D 5 +D ts } 

(cont.) 
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- 0.08696839 



Figure 7f. Expected Rank data of merged 008 strength set {D s +D ts } 

(corn.) 





A 


B 


C 


229 




18.038 


- 0.079584599 


230 


18.064 




- 0.072204977 


231 




18.112 


- 0.064829146 


232 




18.161 


- 0.057456727 


233 


18.183 




- 0.050087343 


234 


18.204 




- 0.04272061 


235 




18.208 


- 0.035356146 


236 




18.227 


- 0.027993566 


237 


18.24 




- 0.020632481 


238 




18.254 


- 0.013272503 


239 




18.286 


- 0.005913238 


240 


18.33 




0.001 445707 


241 


18.351 




0.008804732 


242 


18.405 




0.016164235 


243 


18.421 




0.023524621 


244 


18.46 




0.030886298 


245 




18.464 


0.038249675 


246 


18.49 




0.045615167 


247 


18.49 




0.052983192 


248 


18.51 




0.060354172 


249 


18.553 




0.067728535 


250 


18.566 




0.07510671 


251 


18.586 




0.082489134 


252 


18.599 




0.089876248 


253 


18.602 




0.097268498 


254 


18 . 607 ! 




0.104666334 


255 




18.616 


0.112070215 


256 


18.631 




0.119480602 


257 


18.67 




0.126897966 


258 


18.672 




0.134322781 


259 




18.678 


0.141755531 


260 


18.693 




0.149196704 


261 


18.739 




0.156646797 


262 


18.75 




0.164106316 


263 




18.768 


0.171575771 


264 


18.77 




0.179055685 


265 


18.85 




0.186546586 


266 




18.863 


0.194049014 



Figure 7g. Expected Rank data of merged 008 strength set {D s +D ts } 

(cont.) 
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Figure 7h. Expected Rank data of merged 008 strength set {D s +D ts } 

(corn.) 
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Figure 7i. Expected Rank data of merged 008 strength set {D s +D ts } 

(cont.) 
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Figure 7j. Expected Rank data of merged 008 strength set {D s +D ts } 

(cont.) 
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Figure 8a. Expected Rank data of merged 019 strength set {D s +D ts } 
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Figure 8b. Expected Rank data of merged 019 strength set {D s +D ts } 

(corn.) 
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Figure 8c. Expected Rank data of merged 019 strength set {D s +D ts } 

(cont.) 



77 





E 


F 


G 


115 


14.01 




- 1.127438617 


116 


14.069 




- 1.116643872 


117 




14.085 


- 1.105926941 


118 


14.16 




- 1.095286311 


119 




14.213 


- 1.084720508 


120 


14.27 




- 1.074228096 


121 


14.271 




- 1.063807675 


122 


14.305 




- 1.053457881 


123 




14.328 


- 1.043177384 


124 


14.332 




- 1.032964885 


125 


14.346 




- 1.022819119 


126 


14.353 




- 1.012738851 


127 


14.37 


1 


- 1.002722875 


128 


14.416 




- 0.992770013 


129 


14.434 




- 0.982879117 


130 


14.481 




- 0.973049063 


131 


14.485 




- 0.963278755 


132 


14.528 




- 0.953567122 


133 


14.54 




- 0.943913114 


134 




14.579 


- 0.934315709 


135 




14.586 


- 0.924773904 


136 


14.596 




- 0.915286721 


137 




14.648 


- 0.9058532 


138 


14.649 




- 0.896472405 


139 


1 4.676 




- 0.887143417 


140 


14.678 




- 0.877865339 


141 




14.681 


- 0.868637291 


142 


14.69 




- 0.859458411 


143 




14.697 


- 0.850327856 


144 


14.702 




- 0.841244799 


145 


14.767 




- 0.832208431 


146 


14.825 




- 0.823217957 


147 


14.83 




- 0.814272599 


148 


14.855 




- 0.805371595 


149 


14.929 




- 0.796514195 


150 


14.95 




- 0.787699666 


151 


14.97 




- 0.778927288 


152 


15.003 




- 0.770196353 



Figure 8d. Expected Rank data of merged 019 strength set {D s +D ts } 

(cont.) 
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Figure 8e. Expected Rank data of merged 019 strength set (D s +D ts ] 

(corn.) 
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Figure 8f. Expected Rank data of merged 019 strength set {D s +D ts } 

(cont.) 
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Figure 8g. Expected Rank data of merged 019 strength set {D 5 +D ts } 

(corn.) 
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Figure 8h. Expected Rank data of merged 019 strength set {D s +D ts } 

(cont.) 
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Figure 8i. Expected Rank data of merged 019 strength set {D s +D ts } 

(cont.) 
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Figure 8j. Expected Rank data of merged 019 strength set {D s +D ts } 

(cont.) 
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Figure 8k. Expected Rank data of merged 019 strength set {D s +D ts } 

(cont.) 
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Figure 9. Likelihood surface and contours of 008 strength set {D s } 
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Figure 10. Likelihood surface and contours of 019 strength set {D 5 } 
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Figure 11. 5 % and 1 % confidence intervals for 008 strength set {D s } 
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Figure 12. 5 % and 1 % confidence intervals for 019 strength set {D s } 
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Figure 13. 5 % and 1 % confidence intervals for 008 strength sets {D s } 

and {D ts } 
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Figure 14. 5 % and 1 % confidence intervals for 019 strength sets {D s } 

and {D ts } 
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Figure 15. Expected Rank of 008 strength sets (D s +D ts ) and {D s } 



92 



Strength - s [gm] 



0.999 



(s)*d - aJnjjej p Aijijqeqojd >jun jsq^bqm 



CvJ 



O 



C\J CO LO CO 

i i i i i i 




o 

CO 



CO 



CO 






C\l 



g> CD lo v ^ ro 

^ o d o o o 



(s)d - ©jnijBd jo /aijjqeqojd ©AjieiniunQ 



Figure 16. Expected Rank of 019 strength sets (D s +D ts } and (D s ) 
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Figure 17. 5 % and 1 % confidence intervals for 008 strength sets {D s } 



and {D s + D ts } 
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Figure 18. 5 % and 1 % confidence intervals for 019 strength sets {D s } 

and {D 5 + D t5 } 
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Figure 19. Expected Rank of 008 strength sets (D s +D ts ) and {D^| 5 } 
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Figure 20. Expected Rank of 008 strength sets {D s +D ts } and {D t 2 S } 
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Figure 21. Expected Rank of 008 strength sets {D s +D ts } and {D t 3 S } 
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Figure 22. Expected Rank of 019 strength sets {D s +D ts } and { D t1s ) 
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Figure 23. Expected Rank of 019 strength sets {D s +D ts ) and {D t 2 S } 
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Figure 24. Expected Rank of 019 strength sets {D s +D ts } and {D t 3 5 } 
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Figure 25. Expected Rank of 008 strength sets {D s +D ts ) and (D) 
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Figure 26. Expected Rank of 019 strength sets {D 5 +D ts } and {D} 
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Figure 27. Likelihood contours of 008 strength set {D} 
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Figure 28. Likelihood contours of 019 strength set {D} 
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Figure 29. Likelihood contours of 008 life set {D t] } 
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Figure 30. Likelihood contours of 008 life set {D t2 } 
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Figure 31. Likelihood contours of 008 life set {D t 3 } 
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Figure 32. Strength-Life Model for AS4-008 
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